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Foreword 


This is a study in the general area of multidimensional modeling of Diesel sprays. In this 
work, the vaporization and combustion characteristics of droplets in a Diesel engine are studied 
computationally. It has been determined that individual ‘droplet’ combustion is not relevant in 
Diesel engines. The vaporization rate determines the rate at which vapor fuel becomes available 
for mixing and burning. It also determines whether the liquid impinges on the walls of the 
cylinder. Such impingement is undesirable. The work is carried out from the point of view of 
applications to modeling sprays in Diesel engines. Hence, it focuses not just on vaponzation 
rates of droplets but on determining these drop sizes since the sizes are the most important 
parameter in determining vaporization. The droplet modeling studies are integrated with studies 

of modeling Diesel sprays. 



Statement of the Problem Studied 

Typical liquid injection velocities in a Diesel engine are over 400 m/s. With such high 
velocities, the liquid fuel atomizes into droplets of sizes of 0(1)-0(10) p,m [1-4]. The chamber 
gas temperature and pressure under warm operating conditions are in the range of 900-1100 K 
and 70-120 bars respectively, whereas under cold-start conditions they are 700-900 K and 40-70 
bars respectively. The wide range of drop sizes and thermodynamic conditions in the engine 
chamber makes it necessary to ensure that models for droplet vaporization predict vaporization 
rates with adequate accuracy under all conditions. 

In this work we address the question; Is the current simplified droplet vaporization 
models employed in multidimensional models adequate to capture vaporization rates and drop 
lifetimes? The key unknown variable in modeling the sprays is the droplet size. Hence, this 
work also focuses on modeling the outcomes of droplet collisions that lead to the computed 


sizes. 




Summary of the Most Important Results 


The Computational Model 

Studies are carried out with a detailed one-dimensional droplet vaporization model and a 
simplified droplet vaporization model as employed in multidimensional models [5-7] to obtain 
the results shown here and arrive at the conclusions. 

The One-Dimensional Droplet Vaporization Model 

The governing equations for the model are the conservation equations for mixture mass, 
internal energy and species for both the gas and liquid phases with additional relations for the 
flux terms contained therein [8]. To simplify the problem, the following assumptions are made 
here: 


Al: Droplet vaporization is spherically symmetric. 

A2; All non-radial velocities and fluxes are zero. 

A3: Pick’s and Fourier’s laws for mass and energy diffusion are used. 

A4: Energy terms arising from viscous dissipation, radiation, interdiffusion, and flow 
work are neglected. 

A5: Negligible Soret and Dufour effects. 

A6: Pressure is uniform and constant throughout. 

A7: Thermal, mechanical, and chemical equilibrium at the vapor/liquid interface. 

A8: Lewis and Raoult's assumptions are valid at the vapor/liquid interface. 

With assumptions Al through A6, the conservation equation for mixture mass in both the 
gas and liquid phase is: 


^ + _L i- 

3t 1-2 
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and the conservation equation for mixture internal energy for the gas and liquid phase is: 
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and the conservation of species equation is 
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where p is the mixture density, k is the mixture thermal conductivity, Cy is the constant volume 
specific heat capacity, Dj is mass diffusivity of species “i”, Uj- is the radial velocity, U is 
specific mixture internal energy, and Yj is the species mass fraction. 






In addition, the integral forms of the conservation relations for mixture mass, energy, and 
species are applicable for a control volume which is situated on the vapor/liquid interface. In the 
limit as this control volume becomes infinitesimally thin about the interface the contributing 
terms of the flux components within the plane of the interface become negligible as compared to 
those normal to the interface. Therefore, accounting for the motion of the vapor/liquid interface, 
conservation of mass across the interface yields: 

mis = Pls(uis “ ^s) = ms = Pgs(ugs “ ^s) = mgs (4) 


where pis and pgs are the liquid and gas phase mixture densities, uis and Ugg are the liquid and 
gas phase bulk velocities, and mis and rhgs are the rates of liquid and gas phase mass flux with 
all values being evaluated at the surface and Vs is the interface velocity. 


Conservation of energy across the interface yields: 

Als ~ msAHys + qgs 


(5) 


where qis is the energy flux in the liquid phase, qgs is the energy flux in the gas phase, and 
AHys is the enthalpy of vaporization (i.e., the enthalpy change across the interface) with all 
values being evaluated at the surface. 

And, conservation of species across the interface yields: 

mis(i) = PlsYls(i)(nis + 'Vls(i) “ ’^s) = ms(i) 
mgs(i) = PgsYgs(i) (ngs +^gs(i) ” ^s) =ms(i) 


where the new variables Vis(i) and Vgs(i) are liquid and gas phase diffusion velocities and 
mis(i) and riigsd) are the rates of liquid and gas phase mass flux of species “i” with all values 
being evaluated at the surface. 

From assumptions A7 and A8, the following interface boundary conditions also apply: 

Tis=Ts=Tgs (7) 
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where Tg is the surface temperature, pg is the surface pressure, and Xig(i) and Xgg^j) are the 
mole fractions of species “i” in the liquid and gas phases at the surface, respectively. 

Full account of property variations with respect to pressure, temperature, and 
composition have been included within the model. 

The finite-difference/finite-volume methods of Patankar [9] are used to solve Eqs. (1)- 
(9). The essence of the method involves breaking up the physical domain into a finite set of 
points for which at each point there is a prescribed volume in space, the control volume, 
containing the point. Upon assuming that the dependent variable in question behaves in a 
constant or linear fashion throughout the control volume, the conservation equations are 
integrated over each control volume to obtain algebraic relations in terms of the unknowns. The 
algebraic equations form a system of non-linear, coupled set of equations which are then solved 
in an iterative fashion. In the case of one dimension, the solution may be obtained in one 
application of the well-known tri-diagonal matrix algorithm method. However, due to 
nonlinearities arising from velocity, temperature, pressure, composition dependent properties, 
and fully implicit treatment of boundary conditions, repeated application of the solution process 
is implemented until a converged solution has been obtained. 

SimpliHed Droplet Vaporization Model 

In the simplified droplet vaporization model, the rate of change of mass of the droplet is 
given by the Frossling correlation [6-8]. 



The rate of change of energy of the droplet is given by the equation: 
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In the equations above, rd is the droplet radius, p^ is the liquid density, D is the mass 
diffusivity of the gas, B is the transfer number, Shd is the Sherwood number, C( is the liquid 
specific heat, Td is the droplet temperature, L(Td) is the latent heat of vaporization of the 



liquid, Qd is the rate of heat conduction to the droplet, Yys is the droplet surface mass fraction 
of the fuel, Yy is the mass fraction in the gas and given as Yy = Pu / P. is th® thermal 
conductivity of the gas and Nud is the Nusselt number. 


The following additional correlations are employed to solve these equations: 
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L(Td) is obtained from the Haggenmacher correlation [7] 



where R is the universal gas constant, Wy is the molecular weight of the fuel, P is the vapor 
pressure of the fuel and p is the compressibility factor given by 


P = 




il/2 


where and p^ are critical values of temperature and pressure for the fuel. We have used the 
Atoine equation to obtain the vapor pressure [10] 


B 


logiop = A - 


T + C 



where the constants A, B and C may be obtained for the fuel. 

The properties (pD) , Ag are evaluated at a mean temperature, t, obtained using 
the 2/3 law [6, 11]. 


3 

The simplified droplet vaporization model, described above, is also capable of modeling 
multicomponent droplet vaporization. The multicomponent droplet vaporization model has been 
developed with support from the ARO [12]. 

It has been shown in detailed numerical studies of isolated multicomponent droplets that 
the two important controlling physical processes, which determine the rate at which one 
component vaporizes relative to another, are the volatility of the component as reflected, in part, 
in its latent heat of vaporization and the inter-species liquid diffusivity of the component. The 
greater the volatility and the lower the latent heat of vaporization of the component, the faster it 
vaporizes. The greater the inter-species liquid diffusivity of the component in a multicomponent 
mixture, the faster it diffuses to the surface of the droplet and the faster it vaporizes. Our model 
is built on these essential physics. Consider the physical picture, as shown in Fig. 1(a), of an N- 
component liquid droplet of radius R, with the species mass fractions of .N > in an 

ambient environment of temperature T^, pressure Pg and vapor species Yv^a,i=l.N • There is 

a heart flux of Q into the droplet of which Qy is utilized to vaporize the component and Qi to 
heat the droplet. There are two important characteristic times which then control the 
vaporization process of each species: tv,i is the characteristic time associated with the 
vaporization of the individual species with the more volatile species having a shorter time and 
the less volatile one having a longer time. Td,i is the characteristic time associated with the 

inter-species liquid diffusivity of species I in the multicomponent mixture. The model is 
developed as follows: At any instance in time, as shown in Fig. 1(b), we assume that there are N 
droplets each of a single component species I, radius R and the same ambient environment as the 
multicomponent droplet at the instance in time. The heat flux to each of the single component 
droplet is assumed to be the same as the heat flux to the multicomponent droplet. A vaporization 
rate, rhy j, for each single component droplet and hence each component may then be evaluated 

using the mass and energy balance equations discussed earlier. 

The rate of vaporization of any component in the multicomponent droplet, riiy i, would 

be different from the value computed for this single component droplet, rhy j, since the 

concentration of the species on the surface of the multicomponent droplet will be different from 
1.0. the actual vaporization rate of the species in the multicomponent droplet is modeled as 


itiy,! - ihy j Yj j 


(12) 






where Yj j is the actual component mass fraction on the surface of the multicomponent droplet. 
Yj i will depend on the liquid-phase diffusivity of the species I with the species which has a 
greater diffusivity likely to have a higher concentration on the surface. This is modeled in the 
second of the two essential elements as 


Ym = Yi.i ^ 


Dl,i 




(13) 


where Y\ j is the mass fraction species I in the multicomponent droplet and is the effective 

diffusivity of species I in the mixture. The physical picture represented here is that the 
concentration of the more diffusive species will be greater at the surface than that of the less 
diffusive one and hence more of it will vaporize. In this manner the effects of volatility and 
species diffusion are modeled. 

For the liquid-phase diffusivity of each component, we have used the Wilke-Chang 
equation for the effective diffusivity, Di^,, of each component [10] 


Dij = 7.4 X 10 


,.8 cm2/s 




(14) 


where 
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(15) 


In the equation above, x is the mole fraction of the component in the liquid droplet, M is 
the molar mass, V is the molar volume, <D is the association parameter for the solvent and Ti is 
the viscosity. T|i is obtained from the Grunberg-Nissan equation [10]: 


N 

In Tii = Z Xi In Tii + ^ Z Z Xi Xj Gij (16) 

i=l i j 

where Gy is an interaction parameter which is a function of the components I and j and the 

temperature T. We include temperature dependent properties for the liquid phase density, 
conductivity, specific heat capacity, viscosity and latent heat of formation. 




In the numerical calculation procedure, the equations for droplet mass and energy transfer 
are applied to each of the N droplets, where each droplet is of a single component of the N 
components, at each step and individual vaporization rates are computed. In the case of a droplet 
reaching the critical state, it is assumed to vaporize completely at that point. Once the individual 
vaporization rates are known then the weighted vaporization rates for the components in the 
multicomponent droplet may be computed. An energy and mass balance for the multicomponent 
droplet then leads to updated values of temperature and composition and the procedure is then 
repeated. 

These models were employed in three different geometries which are schematically 
illustrated in Fig. 2. 

The computational grids employed for the detailed droplet computation, the single 
droplet computation using the multidimensional model and the spray computations using the 
multidimensional model are shown in Figs. 2(a), (b) and (c) respectively. Figures 2(a) and (b) 
are schematics. 

Computations have been performed for single droplets of the pure components n- 
hexadecane, n-tetradecane, n-hexane, and benzene vaporizing in quiescent environments initially 
at a uniform temperature, pressure, and composition. Table I lists the critical constants of these 
species. In all studies performed here, the initial droplet temperature was uniform at 300 K and 
initial droplet diameters were 30 pm or 245 pm. 


The Combustion Submodel 

The combustion submodel, that is employed here, includes a model for predicting 
autoignition and one for computing the subsequent combustion of the ignited gas. Our high- 
temperature combustion model follows the evolution of N fuel species and six additional species, 
O 2 , N 2 , CO 2 , H 2 O, H 2 and CO, by solving the follow with i-th species conservation equation 
together with all other flow and spray equations 


8Yi ,, aV; 1 9piVy . 

_^ + Vj -^ =-r—= Wi 

at ^ axj p axj 


(17) 


In the above equation, Y, is the mass fraction of species I, vj is the velocity in the 
direction], xj is the distance in the direction], p is the gas density, Vy is the diffusion velocity 
of species I in the direction ], and wj is the rate of formation of species i. The rate of change of 
species partial density is given by the expression 



'^c 


dPi ^ 
dt 


(18) 



where p* is the local and instantaneous equilibrium value of the i-th species. Xq is a 
characteristic time with which the equilibrium is attained. It is obtained from the following 
expression: 


Tc =max(Tt,T:i) 

where ^ turbulent eddy turnover time and x\ is the characteristic time with which the 
equilibrium is achieved in a laminar flame. %\ is obtained by comparing experimental and 
computed laminar flame speeds. In extending this model to multiple fuel components, the local 
value of the fuel density is taken to be the sum of the partial densities of the different 
components. In computing the local and instantaneous equilibrium species densities, we 
consider the fuel species as a composite one with the formula CxHy 02 where the local values 

of X, y and z may be determined from the local concentrations of the individual fuel species. 

To model autoignition, we add to the 6+N species a new species that represents the 
evolution of the radicals that lead to ignition. The evolution of this additional species is 
controlled again, like that of the others, by convection diffusion and chemical formation. 
However, this species is assumed to be a trace species and hence no mass is associated with it. 
In our previous work, where we considered only one fuel species, the chemical formation was 
due only to that species. In this work, where we have more than one species in the 
multicomponent fuel, the different fuel species will contribute to the formation of this trace 
species and at different rates. In our model, we assume that the activation temperature in the rate 
expression that yields the trace species from each component is inversely proportional to the 
cetane number, CN, of the component [13,14]. The model for the species is then Eq. (17) and 
the following equations 
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gd(l+Cll-(t)|)/(T(CN)i) (22) 

where Dg is the effective diffusivity, i.e., the sum of the laminar and turbulent diffusivities, P is 

the pressure in atm, T is the gas temperature in K, ^ is the equivalence ratio and A, B, C, D are 
constants that are determined by comparisons with measured ignition delay times. 
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Discussion of Results 


Results from the simplified and detailed vaporization models are shown in Figs. 3 and 4 
for single component droplet. The results are for single component hexadecane droplets. Figure 
3 shows results for the computation where the ambient pressure is 10 atm and ambient 
temperatures of 600 K and 900 K are selected. These temperatures encompass the cold start 
operating range. Figure 3 shows that the simplified model reproduces the size history of the 
drops and the lifetime of the drop within 15% for these conditions. 

When the ambient pressure is increased to 50 atm Fig. 4 shows that the maximum 
differences in life history and drop lifetime are still within 15%. This implies that, for these 
conditions, the simplified droplet vaporization model predicts the vaporization characteristics 
with adequate accuracy. Figure 5 shows the drop surface and mass averaged drop temperature as 
a function of time. Provided these temperatures are within 50 K of each other, the simplified 
vaporization model reproduced the detailed results with adequate accuracy. This is certainly the 
case for droplets in Diesel engines under cold start conditions. Figure 6 shows the size history of 
a droplet at an ambient temperature of 1300 K and pressure of 50 atm. In this case the 
differences in lifetime is about 50%. Figure 7 shows the surface temperature and mass averaged 
drop temperature for the droplet under these conditions. Significant differences between the two 
temperatures may be observed. In this case, the simplified vaporization model is not adequate^ 
However, under these conditions corresponding to warm operating conditions in a turbocharged 
engine under heavy load, the liquid vaporization and penetration is controlled by mixing rates 
and not by the rate of droplet vaporization. This has been shown experimentally by [15] and 
computationally [4,16]. 

As part of this work, the conditions under which droplets reach the critical state in a 
Diesel engine have also been investigated. This topic is important, in part, because simplified 
vaporization models are unlikely to reproduce near-critical vaporization with adequate accuracy. 

Figure 8 shows the effect of ambient temperature on droplet size for a constant ambient 
gas pressure of 50 atm. The normalized droplet diameter squared has been plotted as a function 
of the scaled time for a n-hexadecane droplet vaporizing in nitrogen at ambient temperatures of 
600, 700, 800, 1000, and 1300 K. the expected trend of reduced droplet lifetime with increasing 
ambient temperature is reproduced. At ambient temperatures of 800 K and below and after an 
initial transient period, the droplet exhibits the well-known D-squared vaporization regime over a 
significant portion of its lifetime even at this relatively high ambient pressure. However, this 
behavior is not reproduced at higher temperatures. Figure 9 shows the normalized droplet 
diameter squared as a function of the scaled time for the ambient temperature of 1300 K and 
ambient pressure of 50 atm. It may be seen that there no longer appears to be a well-defined D- 
squared vaporization regime throughout the droplet’s history. 

Figure 10 shows the corresponding droplet surface temperatures for the conditions 
considered in Fig. 2. it may be seen that for an ambient temperature of 600 K the droplet surface 
temperature attains the pseudo wet-bulb state for a significant portion of the droplet’s lifetime. 
However, at an ambient temperature of 1300 K the droplet surface temperature is entirely 
transient and, in fact, reaches the critical point near the end of its lifetime. Under these ambient 


conditions of P=50 atm and T=1300 K, the quasi-steady theory would appear to no longer apply 
as seen in both the droplet’s size and surface temperature histories. 

Figure 11 shows the loci of lowest ambient reduced temperature (± 10 K) at a fixed 
reduced pressure which results in the droplet surface reaching the critical temperature of the pure 
component species. Here, the ambient temperatures and pressure have been reduced by the 
critical temperatures and pressures of the liquid species. The critical loci shown for n- 
hexadecane and n-tetradecane are in qualitative agreement with previous suggestions [17] that 
vaporizing pure component droplets will reach the critical state when reduced ambient pressures 
are greater than 2.0 and reduced ambient temperatures are greater than 2.0. In fact, this figure 
shows that for these pure component droplets, when either reduced ambient pressure or reduced 
ambient temperature is greater than 2.0, the value of the other at which the droplet reaches the 
critical state may be less than 2 and its value will decrease as the value of the other increases. 
Figure 11 also shows that the critical loci is fuel dependent. 

The departure of the critical loci for n-hexane and benzene from that of n-hexadecane and 
n-tetradecane may be due to the much greater critical pressures for the former two species, 30.3 
and 48.9 atm, as compared to the latter two species, 13,3 and 14.4 atm, respectively for the sanie 
reduced pressure, the high critical pressure species will be vaporizing in an ambient gas which is 
at a much greater pressure than the low critical pressure species. As a result, the mass diffusion 
rates will be reduced with a corresponding reduction in mass vaporization rates. Since mass 
vaporization is the means by which the droplet carries away the energy flux to its surface, droplet 
heating is enhanced for these high critical pressure species and the droplet will be more likely to 
reach the critical point at a lower reduced ambient temperature. This reasoning would also 
explain why the critical loci for benzene is lower than that of n-hexane. 

Figure 12 shows a comparison of the critical loci for n-hexadecane droplets vaporizing in 
pure nitrogen and carbon dioxide environments. Above a reduced pressure of approximately 2.5, 
there is essentially no change in the critical loci. However, below this pressure the critical loci 
for vaporization in carbon dioxide lies below that of nitrogen. Thus for pressures below 35 atm, 
this suggests that a vaporizing droplet may reach the critical point at a lower ambient 
temperature if the ambient gas contains significant amounts of carbon dioxide. This may again 
be related to differences in the mass diffusion rates of the two species, as explained in the 
previous paragraph, leading to more effective droplet heating in the case of the carbon dioxide 
ambient gas. 

Now, we ask the question: Do multidimensional vaporizing spray models that are in 
widespread use today [6-7] reproduce the trends identified above? In today’s production Diesel 
engines, it may appear unlikely that the compression temperatures reach the values that should 
result in the droplet reaching the critical point. Also, vaporization of the liquid will result in the 
cooling of the cylinder gas. Hence, droplets are likely to reach the critical point only if they are 
within the reaction zones. Experimental results, however, do not support this possibility of drops 
being in the reaction zone [18]. 

To answer the question posed above, next, we present results from computations where 
we employ the multidimensional model [7] to perform single droplet vaporization computations. 



The single droplet is assumed to be at rest in an ambient which is also at rest initially. Figure 
1(b) shows a schematic of the computational domain. The gas void fraction is close to one in the 
computational cell. The liquid species is tetradecane. Figure 13 shows the variation of the 
computed normalized diameter of the drop with time for four cases when the temperature is beld 
constant at 1400 K but the ambient pressures are 20 atm, 32 atm, 48 atm and 68 atm. The initial 

ambient density is kept fixed at 16.67 kg/m^. This is achieved by changing the molecular 
weight of the ambient gas. 1400 K, in the case of tetradecane, corresponds to Tf = 2.02. The 
corresponding Pj for the four cases are 1.39, 2.2, 3.33 and 4.72. On the figures, the abrupt 
transitions in radius results when the droplet reaches a critical state. It may be seen that at Pf = 
1.39, the droplet vaporizes completely without reaching the critical state. However, as the 
pressure is increased, the droplet reaches the critical state at earlier stages in its vaporization. As 
the pressure is raised, the vaporization rate decreases and the self-cooling of the droplet is not 
sufficiently rapid to avoid the temperature reaching the critical value. It is interesting to observe 
that the model predicts that this critical state is achieved at about the same time in all cases. This 
appears to result from a balance between the net rate of heat transfer to the droplet, which 
increases as the pressure is raised, and the greater mass of the droplet whose temperature has to 
be changed at any given time as a result of slower vaporization. These results are consistent with 
those from the detailed computations presented earlier. Figure 14 shows the variation of the 
computed normalized diameter of the drop with time for three cases when the pressure is held 
constant at 48 atm but the initial ambient temperatures are 992 K, 1200 K and 1400 K. The 

initial ambient density is again fixed at 16.67 kg/ m^ in all cases. 48 atm corresponds to Pf — 
3.33 and the temperatures correspond to Tj. = 1.44, 1.74 and 2.02. In this case, the droplet 
reaches the critical state at earlier times and at larger radius as the initial ambient temperature is 
raised. It is seen however that the droplet reaches a critical state even at Tj = 1.44. However, 
this may result because Fj- = 3.33 rather than about 2.0 [17]. Figure 11 shows that at Pf = 3.33, 
Tr for the critical state to be reached is about 1.7. The important conclusion here is that under 
cold start conditions, droplets do not reach a critical state in a Diesel engine. The ambient 
conditions would have to correspond to about twice the critical pressure and about twice the 
critical temperature for the droplet to reach a critical condition. These conditions would 
correspond to warm operating conditions in a turbocharged engine under high load. However, as 
pointed out earlier, under such conditions the vaporization rate is mixing controlled and 
individual droplet vaporization is not relevant in Diesel sprays. 

The single most important parameter that has to be known in order to predict the 
vaporization rates in a Diesel spray is the droplet diameter and this information is not available 
accurately in sprays as it depends on atomization, collisions and drop break-up. As part of the 
work carried out on this contract, a new algorithm that predicts the outcome of drop-drop 
collisions has been developed. The current model includes coalescence and grazing collisions as 
possible outcomes. The new model that has been developed includes bounce, reflexive 
separation, high-pressure effects and drop shattering as additional outcomes. Bounce, reflexive 
separation and high-pressure effects have been modeled to reproduce experimental results of 
[19]. The works of Estrade et al., and Ashigriz and Poo have also been employed to determine 
the different regimes on a Weber number-impact parameter plot. The model that has been 
developed is referred to as the composite collision model. Of these, bounce has been shown to 






have a significant influence on the computed droplet sizes. The work is described in greater 
detail in reference [20]. 

Figure 15 shows the SMD of the drops in a cube (‘box’) where the drops are uniformly 
distributed as a function of time obtained by employing the Brazier-Smith et al. [21] model for 
coalescence efficiency. Results are presented for four resolutions. It may be seen that the results 
do not change significantly for resolutions greater than 40 x 40 x 1. At 100 ms after start of 
computation, the SMD has increased from an initial value of 2 pm to 2.85 pm, an increase of 
greater than 40%. Figure 16 shows the results when the composite coalescence model proposed 
here is employed and compares them with results from the model of Brazier-Smith et al. [21]. 
Results are presented for the same two resolutions where convergence is achieved. At 100 ms 
after start of computation, the SMD has increased to about 2.17 pm with the composite model, 
an increase of less than 10%, whereas it has increased to about 2.85 pm with the Brazier-Smith 
et al. [21] model, an increase of over 40%. Hence, the composite model predicts a noticeably 
lower coalescence efficiency relative to the Brazier-Smith et al. [21] model. The primary 
difference between the two models is the inclusion of sub-model for bounce, which reduces the 
tendency to coalesce. 

We have also examined how the model for shattering collisions affects the computed 
outcomes [22]. Figure 17 shows results with the composite model when shattering collisions are 
not included and compares them to computations when we have allowed shattering to occur for 
collisions with We > 100 in one case and for We > 1000 for another case. The choice of these 
cutoffs for We are somewhat arbitrary and are meant to give insight into how shattering 
collisions might affect computed results. In the computations with shattering, shattering is oitly 
allowed to occur after a rotational or a reflexive separation has been predicted by the composite 
model. Shattering does not occur after a bouncing separating, since the drops to not actually 
touch during bouncing and the disruption of the drops’ surfaces necessary for shattering does not 
occur. As can be seen in Fig. 17, there is negligible difference in results when shattering is 
allowed for We > 1000. This indicates that most of the collisions in this problem are at We 
lower than 1000. When shattering is predicted for We > 100 there is a larger difference initially 
as the drop size reduces by about 10%. The differences arise in the early stages of computation 
when the We numbers are larger and shattering is more dominant. As the initial turbulence in 
the box decays, the We of the collisions becomes lower and the drops stop shattering and 
coalesce. To assess the impact of the coalescence model on computed Diesel sprays, 
computations were carried out in the axisymmetric chamber described in the previous section. 

Summary 

1. A simplified droplet vaporization model that accounts for drop bulk temperature, heat and 
mass transfer effects has been shown to reproduce droplet size history with adequate 
accuracy compared to a detailed droplet vaporization model under cold-start conditions and 
warm operating conditions except under the high pressure (P > 50 atm) and very high 
temperature (T > 1200 K) conditions encountered in a Diesel engine. 

2. It is shown that droplets reach a critical state of vaporization only under the high pressure and 
very high temperature conditions in the engine but certainly not under cold start conditions. 






3. It is shown that under warm operating conditions in a Diesel engine vaporization is 
essentially mixing-controlled and the details of a droplet vaporization model are not relevant. 

4. An improved model for predicting drop sizes in Diesel sprays that result from drop-drop 
collisions has been developed. 




Table 1. Critical Constants of PureRuids 


Fuel 

Tc(K) 

Pc(bar) 

n-hexane 

507.5 

30.3 

n-tetradecane 

693.0 

14.4 

n-hexadecane 

111 

13.3 

Benzene 

562.2 

48.9 



1. PUBLICATIONS SUPPORTED UNDER THIS GRANT 


(a) Peer-Reviewed Journals 

1. V. Iyer, S. Post and J. Abraham, “Is the Liquid Penetration in Diesel Sprays 
Mixing Controlled?” Proc. Combustion Institute, 28: 111 1-1118,2000 

2. J. Abraham and S.D. Givler, “Conditions in Which Vaporizing Fuel Drops Reach 
a Critical State in a Diesel Engine,” SAE Trans., 108: 601-612,2000. 

3. S. Post and J. Abraham, “A Computational Study of the Processes that Affect the 
Steady Liquid Penetration in Full-Cone Diesel Sprays,” Combust. Sci. Tech., 165: 
1 - 20 , 2001 . 

(b) Conference Proceedings 

1. S. Post and J. Abraham, “Modeling the Liquid Penetration in Diesel Sprays,” Proc. 
of the P' Joint Meeting of the U.S. Sections of the Combustion Institute, 224-227, 
Washington, DC, March 1999. 

2. S. Post and J. Abraham, “Challenges in Modeling the Liquid Penetration in Diesel 
Sprays,” Proc. Of the if' Annual Conference on Liquid Atomization and Spray 
Systems, 31-35, Indianapolis, May 1999. 

(c) Manuscripts Submitted 

1. S. Post and J. Abraham, “Modeling the Outcome of Drop-Drop Collisions in 
Sprays,” Int. J. Multiphase Flow (revised), July 2001. 

Scientific Personnel 

Professor John Abraham 

Professor Vinicio Magi (Visiting Professor from Italy) 

Scott Post, Ph.D., 2001 
Shawn Givler 



Bibliography 


1. Bracco, F., “Modeling of Engine Sprays,” SAE Trans. 94,144-167,1985. 

2. Aneja R. and Abraham, J., “How Far Does the Liquid Penetrate in a Diesel Engine: 
Computed Results vs. Measurements?” Combustion Science and Technology, 138, 233-256, 
1998. 

3. Post, S. and Abraham, J., “A Computational Study of the Processes That Affect the Steady 
Liquid Penetration in Full-Cone Diesel Sprays,” Combust. Sci. Tech., 165,1-20,2001. 

4. Iyer, V.A., Magi, V. and Abraham, J., “Exploring Injected Droplet Size Effects on Steady 
Liquid Penetration in a Diesel Spray with a Two-Fluid Model,” Int. J. Heat and Mass 
Transfer, to appear 2001. 

5. Amsden, A.A., O’Rourke, P.J. and Butler, T.D. “KIVA-II: A Computer Program for 
Chemically Reactive Flows with Sprays,” Los Alamos National Laboratory Report No. LA- 
11560-MS, 1989. 

6. Abraham, J. and Givler, S.D., “Conditions in Which Vaporizing Fuel Drops Reach a Critical 
State in a Diesel Engine,” SAE Transactions, 108,601-612,2000. 

7. Magi, V., “REC-2000: A Multidimensional Code for Transient, Two-Phase, Turbulent 
Reacting Flows,” Engine Research Laboratory Report, School of Mechanical Engineering, 
Purdue University. 

8. Williams, F.A., Combustion Theory, Benjamin Cummings Publishing Co., Menlo Park, CA, 
1985. 

9. Patankar, S.V., Numerical Heat Transfer and Fluid Flow, Hemisphere Publishing Co., New 
York, 1980. 

10. Reid, R.C., Prausnitz, J.M. and Poling, B.E., The Properties of Gases and Liquids, McGraw- 
Hill, Inc., 4'” Ed., 1987. 

11. Hubbard, G.L. Denny, V.E. and Mills, A.F., “Droplet Evaporation: Effects of Transients and 
Variable Properties, Int. J. Heat Mass Transfer, 18,1003-1008, 1975. 

12. Abraham, J. and Magi, V., “A Model for Multicomponent Droplet Vaporization in Sprays, 
SAE Transactions, 107, 603-613, 1999. Also in Diesel Combustion Processes, SP-1328, 
117-128, SAE, Warrendale, PA, 1998. 


13. Hardenberg, H.O. and Hase, F.W., “An Empirical Formula for Computing the Pressure Rise 
Delay of a Fuel From Its Cetane Number and From the Relevant Parameters of Direct- 
Injection Diesel Engines,” SAE Paper 790493, 1979. 

14. Olree, R.M. and Lenane, D.L., “Diesel Combustion Cetane Number Effects,” SAE Paper 

840108, 1984. 

15. Siebers, D.L., “Liquid-Phase Fuel Penetration in Diesel Sprays, SAE Paper 980809, 1998. 

16. Iyer, V., Magi, V., and Abraham, J., “Exploring Injected Droplet Size Effects on Steady 
Liquid Penetration in a Diesel Spray with a Two-Fluid Model,” to appear in Int. J. Heat and 
Mass Transfer, 2001. 

17. Givler, S.D. and Abraham, J. “A Review of Supercritical Droplet Vaporization and 
Combustion Studies,” Prog. Energy, Combust. Sci., Vol. 22,1-28,1996. 

18. Dec. J., “A Conceptual Model of DI Diesel Combustion Based on Laser-Sheet Imaging,” 
SAE Paper 970873, 1997. 

19. Qian, J. and Law, C.K., “Regimes of Coalescence and Separation in Droplet Collision,” J. 
Fluid Mech. 331, 59-80,1997. 

20. Post, S.L. and Abraham, J., “Modeling the Outcome of Drop-Drop Collisions in Diesel 
Sprays,” submitted (revised) to Inti. J. Multiphase Flow, July 2001. 

21. Brazier-Smith, P., Jennings, S. and Latham, J., “The Interaction of Falling Rain Drops: 
Coalescence,” Proc. R. Soc. Lond., A 326, 393-408, 1972. 

22. Georjon, T.L. and Reitz, R.D., “Drop-Shattering Collision Model for Multidimensional 
Spray Computations,” Atomization and Sprays 9, 231-254,1999. 





(a) Multicomponent droplet 



(b) Single Component droplets 
Fig. 1. Schematic of model formulation. 




(b) Schematic 



(c) Actual grid 

Fig. 2. Computational domain and grid 
for (a) detailed single droplet, (b) single 
droplet employing multidimensional 
model and (c) multidimensional spray 
computations. 



Fig. 3. Size histories of n-hexadecane 
droplets. 
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Fig. 4. Size histories of n-hexadecane 
droplets. 



Fig. 5. Drop temperature as a function of 
time. 
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Fig. 6. Size history of n-hexadecane 
droplet. 
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Fig. 7. Drop temperature as a function of 
time. 
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Fig. 8. Size histories of n-hexadecane 
droplets vaporizing in nitrogen at 50 
atm. 
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Fig. 9. Size history of n-hexadecane 
droplet vaporizing in nitrogen at 50 atm 
and 1300K. 
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Fig. 10. Surface temperatures of n- 
hexadecane droplets vaporizing in 
nitrogen at 50 atm. 



Pr 


Fig. 11. Critical loci for droplets 
vaporizing in nitrogen environments. 



Fig. 12. Normalized mass vaporization 
rate of a 245 mm C6/C16 (50/50 by 
mass) droplet vaporizing in nitrogen. 



Fig. 13. Variation of drop size with time 
for droplet in ambient with same initial 
temperature of HOOK but different 
pressures. 
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Fig. 14. Variation of drop size with time 
for droplet in ambient with same initial 
pressure of 48 atm but different 
temperatures. 


































































Fig. 15. Computed SMD for 10,000 
parcels uniformly distributed throughout 
a 5 X 5 X 0.5 cm domain at 300K, 1.16 
kg/m^. 1000 drops per parcel. Brazier- 
Smith et al. (1972) coalescence model 
employed. 



Fig. 16. Computed SMD for 10,000 
parcels uniformly distributed throughout 
a 5 X 5 X 0.5 cm domain at 300K, 1.16 
kg/m^. 100 drops per parcel. Brazier- 
Smith et al. (1972) and composite 
models are compared. 



Fig. 17. Computed SMD for 10,000 
parcels uniformly distributed throughout 
a 5 X 5 X 0.5 cm domain at 300K, 1.16 
kg/m^. 100 drops per parcel. Grid 
resolution is 40 X 40 cells. 























